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Cohesive assemblies of filaments are a common structural motif found in diverse contexts, ranging from biological materials 
such as fibrous proteins, to artificial materials such as carbon nanotube ropes and micropatterned filament arrays. In this paper, 
we analyze the complex dependence of cohesive energy on twist, a key structural parameter of both self-assembled and fabricated 
filament bundles. Based on the analysis of simulated ground states of cohesive bundles, we show that the non-linear influence of 
twist derives from two distinct geometric features of twisted bundles: (i) the geometrical frustration of inter-filament packing in 
the bundle cross-section; and (ii) the evolution of the surface geometry of bundles with twist, which dictates the cohesive cost of 
non-contacting filaments at the surface. Packing frustration in the bundle core gives rise to the appearance of a universal sequence 
of topological defects, excess 5-fold disclinations, with increasing twist, while the evolution of filament contact at the surface 
of the bundle generically favors twisted geometries for sufficiently long filaments. Our analysis of both continuum and discrete 
models of filament bundles shows that, even in the absence of external torque or intrinsic chirality, cohesive energy universally 
favors twisted ground states above a critical (length/radius) aspect ratio and below a critical filament stiffness threshold. 
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The superior mechanics afforded by the twisting of ropes and 
fibers has apparently been known to humans for tens of thou- 
sands of years ^ ^. In "staple" yarns and fibers, spun from finite- 
length filaments, twist provides a geometric mechanism of co- 
hesion, transmitting tensile strain along the fibers length to in- 
plane compression of the cross-section, locking the structure 
together via friction^. In continuous strand assemblies, like 
wire ropes and cables, the helical "pre-twist" geometrically 
relaxes the differential stretching/compression under bending, 
conferring a greater flexibility to the structure^. 

Given the robust mechanical properties of twisted filament 
ropes and bundles, it is perhaps unsurprising that Nature in- 
corporates this design motif into a host of biological structures 
at the molecular scale. Helically wound cellulose microfib- 
rils provide mechanical reinforcement for walls of wood cells ^, 
while helically twisted fibrils of extra-cellular protein filaments 
like collagen^ and fibrin^ play crucial mechanical roles in ani- 
mal tissue. Since molecular fibers and bundles of biofilaments 
are not assembled by hand, the spontaneous twist derives from 
the process of assembly itself. Molecular scale chirality of 
biomolecules — the fact that biofilaments are universally heli- 
cal in structure — is broadly implicated as the driving force for 
handed inter-filament packing in twisted bundles^ the fila- 
mentary analogues of the c holest eric and double-twisted tex- 
tures of chiral liquid crystalsEEl. 

While the influence of the twisted geometry of ropes and 



fibers on mechanical properties — high strength under tension 
and extreme compliance to bending — is relatively well de- 
veloped in the textile engineering literature, the question of 
how twist influences inter-filament contact, and therefore the 
cohesive energy in filament bundles, remains poorly under- 
stoodEHm, In this article, we analyze the dependence of pack- 
ing geometry and inter-filament cohesion on the twist of fil- 
ament bundles based on a simple model of inter- filament, at- 
tractive interactions. The dependence of cohesive energy on 
twist is particularly relevant for bundles of molecular (nano- 
)scale filaments, such as carbon nanotubes or filamentous pro- 
teins, where thermodynamic contributions from inter-filament 
attractions may be significant compared to the mechanical costs 
of filament deformation. Even in the simplest model of fil- 
ament cohesion that we analyze here, the geometry of twist 
has a significant influence on both the structure and energetics 
of minimal-energy filament bundles, owing to the surprisingly 
subtle nature of contact between filamentary objects. 

As we show in this study, the complex and non-linear twist- 
dependence of cohesive energy in multi-filament bundles de- 
rives from the interplay from two competing geometrical ef- 
fects: i) geometric frustration of the core packing and ii) sur- 
face energy of non-contacting filaments at the boundary. The 
geometric frustration within the core may be visualized in terms 
of local constraints of packing in twisted filament bundles. A 
measure of these constraints is the kissing number Z, of the 
central filament in a twisted bundle shown in Fig. [T^, which 
counts the maximum number of non-overlapping filaments in 
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Figure 1 (a) Kissing number, Z, for the central (blue) filament, vs rate of twist, Q.d, where d is the filament diameter. Example show the 
structures with the largest Q.d for a given Z. (b) Examples of three cohesive bundles comprised of 145 filaments of a fixed length, L = 40d. 
From left to right, 6 = 0°,3\.2°, and 61.2°. Additional parameters are shown: the filament contour length, L, the bundle radius, R, and the 
twist angle, 6. Coloring is simply used to highlight radial depth of filament position. 



contact with the central filament. The helical shape of the 
finite-diameter filaments surrounding the central filament im- 
plies that contacting neighbors occupy greater than 27r/6 of the 
planar angle surrounding the central (straight) filament. This 
means that at the center of the bundle, twist obstructs the six- 
fold hexagonal close-packing that is possible for straight and 
parallel fibers. Therefore Z < 6 for any finite degree of twist, 
measured by a non-zero rate of twist, Q., with 271 /Q. defined as 
the helical pitch of the bundle. The impossibility of perfect six- 
fold lattice packing at the core has recently been shown^^ to 
derive from an exact mapping of the contact geometry of fila- 
ments in a twisted bundle on to the problem of packing of discs 
on a spherically-curved surface, whose Gaussian curvature is 
proportional to Q?. This mapping implies that, like the ordering 
of particles on spherically-curved surfaces 1) the geometric 
frustration of filament spacing introduces inter-filament strains 
whose cost grows with bundle twist, and 2) in sufficiently large 
and twisted bundles, an excess of certain topological defects, 
5 -fold disclinations, a re nec essary components of the minimal- 
energy lattice packing^ Hence, a quantitive analysis of the 
cohesive energy of twisted filament bundles requires careful ac- 
counting of the number and organization lattice defects, whose 
presence alters the inter-filament spacing throughout the bun- 
dle. 

While the geometric frustration of inter-filament packing is 
largely localized near the interior core of twisted bundles, the 
energy cost of non-contacting filaments is distributed at the 
outer surface of the bundle, whose geometry is also strongly 
altered by twist. In untwisted bundles of uniform length fil- 
aments, the surface energy derives only from the cylindrical 



sides of the bundle, as the ends expose no length of non- 
contacting filaments. As a bundle twists, the constant length 
and helical shape of filaments implies that the vertical height of 
the bundle shortens with radial distance from the bundle cen- 
ter, resulting in a "tapering" of the bundle profile that increases 
with Q,. As the occupied volume within the bundle is largely 
conserved, the tapered bundle must also expand radially, lead- 
ing to a complex progression of surface profiles with increas- 
ing twist, as shown in Fig. [TJ). Due to this complex boundary 
shape, determining whether the surface energy of the bundle 
increases or decreases with twist becomes a delicate affair. The 
ends of twisted bundles expose non-contacting filament length, 
whereas the decreased area of the cylindrical sides reduce this 
exposure. The overall balance of surface energy depends criti- 
cally on just three length parameters of the bundle: the filament 
contour length, L, the (untwisted) bundle radius, Rq, and the 
twist angle, 6, which is the angle between the outermost fila- 
ments with the vertical twist axis. 

In this article, we analyze the thermodynamic consequences 
of both packing frustration and surface geometry in twisted 
bundles based on a combination of numerical and analytical 
studies of a simple model of short-ranged cohesive interactions 
between finite-length filaments. In ref. we introduced and 
explored a simple model of cohesive interactions in A/^-filament 
twisted bundles to study the dependence of the minimal-energy 
filament packing on twist. Here, we extend the analysis of 
this model, by studying the evolution of the cohesive energy in 
with twist, accounting specifically for changes associated with 
filaments at the sides and ends of finite-length filament bun- 
dles. We show that this twist-dependence can be decomposed 



2 



into two contributions: 1) a bulk, elastic contribution that ac- 
counts for the effects of frustration of inter-filament spacing in 
the core, and the subsequent lattice defects which relax frus- 
tration; and 2) a surface contribution, accounting for the total 
length of non-contacting filaments in the bundle. We develop a 
continuum theory of the surface energy, and show that the de- 
crease of surface area at the bundle sides by twist dominates the 
cost of increased filament exposure at the ends, provided a large 
enough aspect ratio L/R^. Importantly, this analysis shows that 
the relaxation of surface energy generically dominates the me- 
chanical cost of twist-induced filament bending, provided that 
inter-filament cohesion (per unit contact length) exceeds a uni- 
versal critical value of 70 > l.297B/Rod, where B is the bend 
modulus of filaments and d is the filament diameter. Based on 
these results and the results of our numerical model, we are led 
to the surprising conclusion that — even in the absence of intrin- 
sically chiral interactions driving twist — sufficiently long and 
flexible filaments generically favor a twisted geometry, pro- 
vided that bundle packing adopts the distribution of lattice de- 
fects required to mitigate the frustration at the core. 

The remainder of this article is organized as follows. In sec- 
tion [2] we introduce the discrete model cohesive interactions in 
A^-filament, twisted bundles. In section |3] we analyze the struc- 
ture and energetics of ground-state packings of bundles in the 
limit of infinite length filaments. In section [4j we develop a 
continuum analysis of the surface energy in finite-length bun- 
dles, which we compare to the infinite-length results. Addi- 
tionally, we analyze the thermodynamics of surface energy in 
finite-length twisted bundles in combination with the mechani- 
cal costs of filament bending. In section [5] we consider the en- 
ergetics of our discrete model for the case of finite-length bun- 
dles, analyzing the thermodynamics while accounting for both 
core and surface packing of filaments, and determine the con- 
ditions for which minimal-energy bundles have non-zero twist. 
We conclude in section [6] with a discussion of the importance of 
filament packings on the cohesive energy, as well as an analysis 
of the predictions of our model in the context of a range of real 
cohesive filament systems. 



2 Discrete model of filament cohesion in 
twisted bundles 

In this section we introduce a model of A/^-filament twisted bun- 
dles possessing short-range cohesive interactions. First, we de- 
rive the form of the pair- wise contact interactions between fila- 
ments in terms of the local geometry of inter-filament contact. 
We then analyze the specific geometric effects on contact of 
filaments within the bounds of a twisted bundle geometry. 



2.1 Filament contact and cohesive interactions 

Here, we summarize the relationship between pair-wise inter- 
actions between portions of filaments (arc-length elements) and 
the cohesive energy expressed in terms of the local geometry 
of inter-filament contact. In this study, we consider bundles of 
homogenous and mutually- attractive filaments, where the in- 
teractions and mechanics can be described purely in terms of 
the shape of the filament center line, X; (^/), which describes 
the position of filament / at arc-position Si along backbone. 
We assume that the filament pairs interact via the summation 
of short-ranged, pairwise interactions between arc-length ele- 
ments, such that the interaction between filaments / and j has 
the form, 

Eij = J dsij dsjy{\Xi{si)-Xj{sj)\), (1) 

where V (r) is an isotropic, finite-range potential describing in- 
teractions between length elements separated by a distance r. 

While a given length element at Si on the filament / interacts 
with the entire length on j, the finite range of segment inter- 
actions generically implies that interactions of filament / at Si 
are dominated by region of filament j closest to ^i{si), which 
we call the contact region. Sufficiently, far from the filament 
ends, we may describe the contact geometry of Si with filament 
j straightforwardly, in terms of s'j{si), a function that maps the 
arc-position on / to the position on j closest to X; (^;), which 
we call the point of contact. Notably, this allows for a generic 
expansion of the shape of the contacting filament around the 
point of contact, in terms of the local geometry of j and dis- 
tance along j from this point, 5s j = sj — s'j{si), 

Xjisj) = Xj{s*) + SsjTj + (5.,)2|n,- + O [{5s j)'] , (2) 

where T^-, Nj and Kj, are the tangent, normal and curvature 
of filament j at ^^(^j) From this expression we find the 
interfilament distance, A = ^j{sj) — ^i{si), from 

I A(5.,-) p = I A,-,f + (5.,-)' (1 + ^Aij • Nj) + O [{8s jf] , (3) 

where A^j = Xj(^p — X;(^;) is the distance of closest contact 
from Si to filament 7, such that A^^ • = 0. These parameters 
are shown schematically in Fig.[2j 

Assuming that V (r) is sufficiently short-ranged compared to 
filament length and curvature, we may use eqn ([3]) to derive the 
cohesive energy contribution of the length element at Si, 

dEij = dstJ^y{5sj)V{\A{5sj)\) = -^^M^=.dsi, 

(4) 

where 

riAij) = j^^du V (^4+^) , (5) 
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is a local contact potential, describing the cohesive energy per 
unit length of filament / in contact with filament j. This po- 
tential is a function of the local distance of closest contact be- 
tween the filaments. The derivation of this local inter-filament 
cohesive energy is fully general for any short-ranged poten- 
tial. In the remainder of the article we study the case of a 
Lennard- Jones interaction between length elements, V(A) = 
£[(cr/A)^^ -2(cr/A)^], for which the contact potential be- 
comes, 
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which like the Lennard- Jones, has stiff repulsion at short-range 
and soft attraction at long range. This potential is characterized 
by an attractive well of depth 70 = 1.6862ecr, at a preferred 
spacing d = 0.947 Icr, which we denote as the filament diame- 
ter. 

Eqn ^ describes how short-ranged, pair-wise interactions 
between all length elements of a filament pair can be formu- 
lated in terms of functions of the local contact geometry (e.g. 
Aij, Kj, and Nj) integrated over the contacting length of a single 
filament (in this case, filament /). In the following section, we 
derive the contact geometry for filament pairs in twisted bun- 
dles. 



2.2 Filament contact geometry in twisted bun- 
dles 

In this section, we derive the geometry of inter-filament con- 
tact in bundles twisted around the central z axis at a uniform 
rate of rotation Q,, which has been presented in ref. We 
consider bundles whose packing is homogeneous along their 
length, such that cross-sectional packing in a plane perpendicu- 
lar to z differs only by the rigid rotation about the z axis. Defin- 
ing the position of filament / at a common plane z = (arc- 
position Si = 0), by the polar coordinates {pi^(j)i), the shape of 
the filament backbone follows the helical curve, 

^i{zi) = Pi cos + Q,zi)x + Pi sin + Q.Zi)y^ ZiZ, (7) 

where it is convenient to express position in terms of vertical 
height Zi = Si cos 6 (pi), where 6{p) = arctan(^lp)is the helical 
angle of filament / with respect to z, shown in Fig.|2p[^ The unit 
normal and curvature are easily calculated from the second- 
derivative of Xj with respect to arc-length. 



Nj = — cos ((^j + Q.Zj)x - 



sin [(j)j + Q.Zj)y. 



(8) 
(9) 



^Note, that 6{p) refers to the local tilt angle of filaments at a radius p in the 
bundle, where as, in our notation, when tilt angle appears without an explicit 
radius, as 0, it refers to helical angle at the outer radius of the bundle, which 
we call the twist angle of the bundle. 




Figure 2 (a) One and a half pitch lengths of a contacting helical 
filament pair, / and j. The helical radius and helical angle for filament 
j are pj and Oj{p) respectively, (b) A view of a region of contact, 
showing the filament parameters defined in eqns |2]) and ([3|. 



For a pair of filaments, / and 7, it is convenient to express the 
separation between curves in terms of a vertical offset, Zij = 
Zi — Zj, and the angular separation in the plane, = 0^ — 0^, 



^^{zij)=pf^pj-2piPj cos {(t)ij^azij) -^zl 



(10) 



The distance of closest contact between / and j is determined by 
the minimization of A^{zij), with respect to the vertical offset 
between contacting points, to find the height separation at the 
distance of closest contact, zjj. From dA^{zij) /dzij = we find 
a transcendental equation satisfied by z^j 



-Q}piPjsin{^ij^Q.zl 



(11) 



To make practical use of this condition, we examine the so- 
lutions in the limiting cases where Q?pipj is either small or 
asymptotically large. The former case corresponds to either 
small twist, or filament positions sufficiently close to the center 
of the bundle, in which case the Taylor series expansion of eqn 
( pT) yields. 



lim ^Iz*/ = 



Q?'pipj sin (\)ij 
\ ^Q? Pip j cos ^i^ 



-0[{Q}piPj)% (12) 



which represents a modest tipping out-of-plane of the distance 
of closest contact between azimuthally- separated filaments. In 
the opposite limit, filaments are far from the core of the bundle 
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in comparison to their helical pitch. Since z *y is always strictly 
less than the pitch, the left-hand side of eqn |TT] ) is never larger 
than In in magnitude, hence, in the limit 0?pipj 1 , the solu- 
tion for Zij must satisfy sin(0/j oc {Q^pipj) so that in the 
asymptotic limit of large helical angles, we have 



based on the approximation of A/^ = A(z*^) described previ- 



lim Q.zlj = -^ij, 



(13) 



which shows that the inter-filament distance of closest contact 
is nearly vertical far from the core of bundle. 

These simple results can be incorporated into an approxi- 
mate formula for vertical contact separation that interpolates 
between the two asymptotic limits. 



- arctan 



(t 



Q.^ Pip j sin 
^ Q? pip j cos (i^i 



(14) 



We employ this form of zlj to approximate distance of closest 
contact in our numerical studies via the relation A^^ = A(Z;*^). 
Notably, this approximation of A^^ becomes poor only for fil- 
ament pairs distant from the center (large p) and on opposing 
sides of the bundle {^ij ^ tt), thus providing an accurate de- 
scription for filament interactions in twisted bundles incorpo- 
rating sufficiently short-ranged potentials {\0^d\ <C 1). 



3 Discrete model simulations: Core 
packing and energetics of twisted- 
bundle ground states 

In this section, we study the ground-state packings of twisted 
filament bundles by numerical simulation of our discrete 
model. Here, we expand upon the study of the evolution of 
filament packing with bundle twist presented in ref. corre- 
lating non-linear structural changes of the ground state with 
the complex twist dependence of the cohesive energy. As 
we find below, the appearance of lattice defects in the cross- 
sectional packing underlies certain universal features of the 
twist-dependence of energy. We first review our simulation 
method, and then present the results of simulations for defect 
structures of ground-state packings. We then proceed to ana- 
lyze the generic form of the twist-dependent cohesive energy 
in our discrete model and its connection to the structural evolu- 
tion of optimal packing. 

3.1 Numerical simulations of twisted bundle 
ground states 

To study the twist-dependence of the ground- state packing and 
cohesive energy of filament bundles, we numerically minimize 
the sum of pair- wise interactions described by eqns (|4]) and ([6]), 



ously in eqn ( 14 ). Since the interfilament contact geometry is 
constant along the bundle height, the cohesive energy for the / 
and 7 pair becomes. 



Eij = 



(15) 



where Lij is defined as the length of filament / in contact with 
filament j. In this section, we consider only the case of infinite 
length, where filament contact is maintain along the full con- 
tour Lij = L, provided that we take / to be the "outer" filament 
(pi > pj). In this way we may neglect any explicit change of 
filament contact length at the ends of the bundle in the L ^ oo 
limit. We numerically optimize the total cohesive energy for 
bundles of fixed N and Q.. Our approach, described in detail 
in ref. generates 500 — 3000 random initial configurations of 
in-plane filament positions, {X;(^/ = 0)}, and seeks a minimal- 
energy configuration of in-plane filament positions using the 
method of steepest descent. From this, we have found the min- 
imum energy packings for bundles with N between 16 and 196, 
and a multitude of twist angles between 0° and ~ 80° . 

3.2 Defects in filament bundles 

The numerical simulations described above produce the 
ground-state positions and energies of filaments in the cross 
section of a bundle, given the input of N and Q. These resul- 
tant packings are analyzed in terms of the network of nearest- 
neighbor contacts. The bond network is formally defined 
through a specialized Delaunay triangulation that takes into ac- 
count the discrepancy in the separation between filament po- 
sitions in the xy plane, and the true distance of closest con- 
tact between filament pairs The triangulation of positions 
produces a unique network of non-overlapping bonds, whose 
"nearest neighbor" connections serve as a measure of the local 
packing geometry. While in the case of an untwisted bundle 
the network universally describes the 6-fold, hexagonal pack- 
ing, for sufficiently large twist, uniform 6-fold coordination is 
not maintained and the packing becomes interrupted by topo- 
logical defects. 

A primary type topological defect in 2D crystalline pack- 
ings is a disclination, which describes the breakdown of long- 
range ^-fold orientational symmetry in the lattice at a singular 
point In a hexagonal crystal, low-energy disclinations are 
typically of two types: 5-fold and 7-fold. At the core of a 5- 
fold disclination, is a lattice site (corresponding to a vertex in 



^As described in ref. to determine the contact network, in-plane filament 
positions are transformed via a conformal mapping in the plane. This mapping 
has the property that (infinitesimally) small in-plane separations in the mapped 
coordinates along the radial and azimuthal directions correspond to the same 
interfilament distance. The nearest-neighbor bond network is computed by a 
Delaunay triangulation of these mapped coordinates that correctly account (lo- 
cally) for the distances of closest contact rather than the in-plane distances. 
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the bond network) with five nearest neighbors, one fewer than 
the six neighbors of the perfect lattice. A discHnation may be 
introduce in a hexagonal lattice by cutting out and removing a 
60° wedge, and stretching it back to close the gap. While the 
large elastic strains generated throughout the crystal generically 
make even single disclinations prohibitively expensive in most 
types of (planar) 2D crystals, our simulations reveal that discli- 
nations are necessary components the minimal-energy pack- 
ings of twisted filament bundles. 

Following the topological characterization of disclinations in 
crystalline solids we define the total topological charge of 
disclinations as 

Q = Y,{6-n)V{n), (16) 

n 

where ^ is a coordination number of nearest neighbors belong- 
ing to a filament, and is the number of internal filaments 
that possess n nearest neighbors. Non-internal, or boundary fil- 
aments, are distinguished by having at least one neighbor bond 
on the outer edge of the bond network. 

Throughout this aritlcle, we consider the structure and ener- 
getics of three discrete bundle sizes in details: small (A^ = 34); 
intermediate (A/^ = 82); and large (N = 184). In Fig.|3]we show 
the evolution of ground- state packings these three bundle sizes. 
The packing of each of these bundles exhibits a common trend 
with increasing twist. Packings evolve from defect-free 2 = 
at zero twist, to an increasing value of topological charge — 
characterized a universal excess 5-fold defects with increasing 
twist. As the intermediate and large bundle packings illustrate, 
while the ground-state packing may include 7-fold disclina- 
tions, negatively-charged defects are always sufficiently out- 
numbered such that the net charge increases with twist. 

By analyzing the topological charge, we construct a phase 
diagram bundle ground states in terms of twist angle, 6, and 
number of filaments, N, shown if Fig. [4^. Notwithstanding 
some obvious limitation of small bundle sizes to reach large 
values of Q, it becomes clear that the optimal value of topo- 
logical charge is predominantly determined by 6, and largely 
independent of, N. Q reaches a maximum of six at > 70°. 
Further simulations show that Q does not increased above 6, 
even for twist angles of up to 87°. We discuss the geometric 
origin of this universal dependence of 2 on in the next sec- 
tion. 

While the net charge is largely fixed at a given twist angle, 
there are many ways to achieve a particular value of Q. Specif- 
ically, 5 and 7-fold disclinations that appear in "neutral" pairs 
do not adjust this value. A commonly observed example oc- 
curs for 2 = 1 bundles, which have n 5 -fold disclinations and 
(^ — 1) 7-fold disclinations. For large bundles this feature be- 
comes important, as shown in the 2nd column of Fig.[3](bundles 
Ila, lib, and IIC). While these three examples have the same 
twist angles and hence maintain Q = I, the total number of 
disclinations, Noisdination^ increases from 1 to 3 to 5 for small. 



intermediate and large bundles, respectively. 

A second phase portrait of the number of disclinations, 
^Disclinatiom P^r charge, 2, is show in Fig. Unlike the 
net disclination charge, the optimal value of Noisdination/Q in 
ground-state bundles varies with both N, as well as 0. In Fig. 
|4]3 we have roughly delineated the 2 7^ behavior into four 
regions. In the small-A^ red region, Noisdination = 2, and bun- 
dles contain primarily only 5 -fold disclinations, without excess 
5-7 pairs. For larger filament number, shown as a light blue 
region, Noisdination = 32- Ground states within this region con- 
tain structures such as lib, where there are one 7-fold and two 
5 -fold disclinations per topological charge. At even larger bun- 
dle sizes is the pink region that continues this trend, i.e. there 
are two 7-fold and three 5 -fold disclinations per topological 
charge, as seen in lie. This trend continues with increasing 
A^, reaching Noisdination = 7 (as seen in Fig.|4|:) and higher for 
simulations not shown. A final region is loosely defined at the 
upper limits of twists, where packings identified in our ground- 
state search belong to a more complex taxonomy, possessing in 
general, non-integer values of Noisdinatwn/Q- This is achieved 
using solely 5 and 7-fold disclinations, like in Xb, Xc, and Xd, 
or with higher-charge disclinations like the 8-fold disclination 
shown in Xa (Fig.|4]l). 

3.3 Mapping between twisted bundles and 
curved surfaces 

In the previous section, numerical simulations of twisted bun- 
dles revealed a complex structural evolution of ground- state fil- 
ament packings with progressively increasing twist, marked by 
an increasing number of lattice defects in the cross-sectional 
order. The complex spectrum of defects is characterized by a 
universal dependence of the topological charge of the packing 
on the twist angle. In the next section, we show that transitions 
in the topological charge of the packing correspond to sharp, 
singular features in the twist-dependence of cohesive energy. 
As the twist-dependence of topological charge underlies a key 
feature of bundle energetics, in this section we briefly review 
the connection between twisted bundle packing and packing on 
spherically-curved surfaces, first described in the ref.^. 

The universal twist dependence of Q in ground- state bundles 
derives from a formal mapping of the inter-filament distance 
and the geodesic distances between points on a dome-like sur- 
face that encodes the metric properties of bundles. The de- 
tails of this mapping can be found in Appendix [A[ To summa- 
rize, the positions of filament centers in the cross section of a 
twisted bundle can be mapped exactly to positions on what we 
will refer to as the bundle -equivalent surface. This surface is 
spherical at the top — corresponding to the small 6 center of the 
bundle — and tapers to a cylindrical geometry along its height — 
corresponding to the large 6 exterior of the bundle. Starting at 
the center of a bundle and traveling along an outward radial 
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Figure 3 Example cross sections of for three example bundle sizes and seven example twist angles, 0, showing two trends: 1) the topological 
charge, Q, increases with twist angle, and 2) the total number of defects, NDiscimation^ increases with bundle size. 5 and 7-fold disclinations are 
colored red and blue respectively. The bottom left two cross sections show their triangulated nearest neighbor bond network. The roman 
numeral labels match specific bundles to their locations within Fig.|4] 



path, is equivalent to starting at the top of the dome and travel- 
ing outward along the surface. An example of a bundle packing 
mapped in this way is shown in Fig. [5] 

The mapping between the interfilament distance in twisted 
bundles and the geodesic distance on the dome-like surface, 
implies that the packing of the bundles can be understood in 
terms of the packing of point-like objects on the a curved sur- 
face, for which classical theorems of differential geometry re- 
late the geometry of a surface to the topology of the nearest- 
neighbor bond network A finite-diameter bundle corre- 
sponds to a finite "crystalline" patch of particles or discs ex- 
tending over a portion of the dome-like surface, centered on 
the high-curvature central pole. In Appendix [A| we derive a 
simple generalization of the Euler-Poincare theorem for a tri- 
angular bond-network on a curved 2D surface that relates the 
topological charge Q of the packing to the integrated Gaussian 
curvature over the surface domain covered by the patch. Based 
on the simple assumption that large deviations from equilateral 
packing at the boundary are energetically expensive, and there- 
fore, unlikely in ground-state packings, we exploit this theorem 
to derive an expression for Qi^, the ideal value of topologi- 
cal charge of the interior packing, which requires no distortion 
from equilateral bond-order at the patch edge (and equivalently, 
the bundle surface). Assuming an axisymmetric shape for the 
boundary of the packing, evaluating the integrated curvature on 
the dual surface gives the ideal charge purely in terms of twist 



angle, 

aw = 6(l-cos^e). (17) 

Importantly, the -dependence of Qid encodes the increase of 
integrated Gaussian curvature as the patch size grows large 
compared to the curvature radius of the dual surface (propor- 
tional to the pitch). As the surface domain grows to cover 
a larger portion of the curved surface, the preferred topologi- 
cal charge becomes non-zero. Though the presence of the free 
boundary of the bundle allows the topological charge to adjust 
based on purely energetic considerations, the positive curvature 
of the bundle-equivalent surface suggests a connection between 
the optimal packings of highly twisted bundles and the better- 
known packings of particles (or discs) on closed, spherical sur- 
faces, studied in the context of the generalized Thomson prob- 
lem——. In the language of the continuum theory of curved 2D 
crystals, excess 5 -fold disclinations screen the elastic stresses 
generated by geometric frustration, such that increasing cur- 
vature (or twist in the case of bundles) requires an increasing 
number of "neutralizing" disclinations. 

Though Qi^ increases continuously from at = 0° to 6 at 
= 90°, the actual topological charge of the packing may only 
take on integer values, the simple assumption that the integer 
value of Q in ground-state packings is determined by the clos- 
est integer value to Qid{0) is remarkably consistent with our 
numerical simulations of twisted bundles. In Eig.[6]we plot the 
Q for bundle simulations vs. twist angle and compare this to 
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Figure 4 (a) Phase diagram of net disclination charge, Q, with 6 and N. (b) Phase diagram of NjyiscUnationsI with twist angle 0, and number 
of filaments, N. Black lines are shown to roughly delineate regions of qualitatively distinct ground-state packing. The roman numeral labels 
correspond to bundle cross sections seen in Fig.[3]and part (d). (c) A zoomed-in view of a seven disclination long Q=\ scar, with bond 
network shown, for a bundle with N =\96 and Q = 25.7°. In (d) we show five examples of bundles that exist in the variable-size defect 
region. 5, 7, and 8-fold disclinations are colored red, blue, and purple respectively. 
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Figure 6 The ideal charge, g/^, from eqn ([TtJ (soUd black line) vs 
twist angle 6. The data points are the Q values for an accumulation of 
5000 simulation results for the various twist angles and bundle sizes 
introduced in section [TT] The number of filaments range from 
N = \6 (red data points) ioN = 196 (blue data points). 



Figure 5 (a) A bundle with 70 filaments and 6 = 57.7°. (b) Top and 
side views of the corresponding disc packing on the 
bundle-equivalent surface, which shares the packing geometry and 
topology with the twisted bundle. 5 and 7 -fold coordinated elements 
are colored red and blue respectively. 



the continuously increasing value of Qid{0). While the agree- 
ment is between Q and Qtd is imperfect, this above argument 
highlights the fundamentally geometric nature of packing frus- 
tration in twisted bundles, and importantly, provides a natural 
explanation for the observation that optimal values of Q are in- 
dependent of N, determined only by 6, which controls the value 
of the integrated curvature on the bundle-equivalent surface. 

This mapping of filament packing in bundles to a spherically- 
curved surface not only explains the increase of Q with twist, 
but it also sheds light on a basic tendency of increasing the 
number of disclinations per topological charge with bundle 
size, as observed in our simulated ground states. In struc- 
tures where Noisdination = Q-\-2n, a defect "complex" can be 
described as a charged scar, usually observed as chain of al- 
ternating 5 and 7-fold disclinations. Importantly, the transition 
from compact (point-like) to extended (string-like) disclination 
defects has been predicted and studied in great detail in the con- 
text of crystalline order on spheres . Adapting the scaling 
arguments for stability of dislocations (neutral 5-7 pairs) in the 
vicinity of a scar on a sphere we expect the number of 5-7 
pairs decorating a scar in a bundle to grow as ^ Ro/d oc ^/N 
at fixed value of twist, a trend that is in reasonable agreement 
with our simulations. 



3.4 Twist-dependence of cohesive energy: 
Infinite-length Hmit 

Having analyzed the twist-dependence of the structure of 
minimal-energy bundles in terms of the topological charge of 
the packing, we now consider the twist-dependence of cohesive 
energy, with the aim of discerning the influence of the universal 
evolution of Q on energy. We find that in the laigc-N limit, bun- 
dle energetics converge to a common behavior. Here, we focus 
on thermodynamic behavior of small, intermediate and large 
bundles in detail. Fig. [7] shows the change in mean cohesive 
energy per filament length for these three cases. We observe a 
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Figure 7 Mean filament cohesive energy per unit length vs twist 
angle, 6, for three selected bundle sizes: small N = 34 (red), medium 
N = S2 (blue), and large N =1S4 (green). The roughness is a 
consequence of sudden rearrangements in the packings to 
accommodate the twist-induced geometric frustration. In this and all 
following plots, the energy change is defined relative to the untwisted 
case. 
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common behavior for each of the bundle sizes. Twisting ini- 
tially increases the energy of untwisted bundles, until reaching 
a rough plateau region at intermediate twist. However, further 
twist lowers the cohesive energy, ultimately driving it below the 
energy of the untwisted state. 

In the next section, we carefully analyze the surface geom- 
etry of twisted bundles to show that the tendency to decrease 
the cohesive energy with twist is driven by a decrease of non- 
contacting filament lengths at the boundary of long bundles. 
In this section, we show that removing the effects of changes 
in filaments at the boundary of the bundle reveals a universal 
dependence of the core packing on twist. To perform this anal- 
ysis, we subtract the surface energy of the bundle, which is 



defined formally in section 4.2 from the total energy of the 
discrete model to define the bulk cohesion energy E^uik- Fig- 
[8] shows the change in mean E^i^ik per filament length for large 
bundles in the range ofN= 166 — 193, revealing a common 
increase of energy relative to the untwisted state. For small 




Figure 8 Bulk cohesive energy vs. twist angle for ten values of N 
from 166 to 193, showing convergence a of internal energy for large 
N. The discrete series of minima correspond to discrete increases in 
Q. Background colors correspond to the trends in Q presented in Fig. 
|4] Inset shows three bundle sizes: small N = 34 (red), medium 
= 82 (blue), and large = 184 (green). 

twist, the core energy increases smoothly with 6, due to the 
increasing frustration of inter-filament spacing in defect free 
bundles. The small- dependence of large-A/^ bundles is con- 
sistent with results of elasticity theory calculations which show 
that Eiji^ik ^ 0^ in defect- free bundles^^^^. Also consistent with 
elasticity theory results, is the appearance of cusps in Ehuik at 
transitions in topological charge, such as the transition from 
2 = to 2 = +1 at 27°. Just beyond the transition, en- 
ergy decreases with twist, highlighting the key ability of excess 
disclinations to mitigate the twist-induced frustration, taming 
the rapid growth in bundle energy. The five local minima mark 
values of twist that are locally stable due the optimal screening 
provided by discrete values of 2 = 1 , 2, 3, 4 and 5. 

We consider the convergence to the large-A/^ core-energy in 



the inset of Fig. [8] which shows the twist dependence of Eij^ik 
(per filament) for the small, intermediate and large bundles. 
Significantly, at a given value of twist, the bulk per filament 
cost of twist decreases as N increases. Underlying the N- 
dependence for small bundles are two distinct types of ener- 
gies counted in ^^^^/yt- The first contribution is the far-field 
elastic cost of geometric frustration and excess disclinations 
in twis ted bun dles. Based on the continuum theory analysis 
of refs.lSElES^ these -dependent costs are proportional N and 
determined only by Q, consistent with the universal large-A^ 
behavior of E^uik- The second type of contribution may be 
attributed to the "core energy" of the excess disclinations ap- 
pearing in the packing above the critical twist. Core energies 
represent corrections to the continuum theory behavior due to, 
for example, the fewer cohesive contacts at the core of 5 -fold 
disclinations, as well as the large- strain deviations in the neigh- 
borhood of defects. These costs grow in proportion to the total 
number of defects in the bundle, and hence amount to a larger 
proportion of E^uik for smaller bundles. As the total number of 
defects does not grow faster than A^, the defect cores contribute 
a negligible amount to the energy density Eij^ik/N relative to 
the far-field elastic costs of disclinations and twist as A/^ ^ oo. 

To summarize, this analysis of the bulk cohesive energy of 
ground-state packings of our discrete filament model reveals 
two key influences of twist on the cohesive energy of bun- 
dles. First, twist frustrates the uniform inter-filament packing 
allowed in straight bundles, leading to a necessary increase in 
energy with twist angle. Second, we find that increasing twist 
triggers the stability of excess 5 -fold disclinations in the cross- 
sectional packing that mitigate the growth of bundle energy 
associated with twist-induced frustration of nearest-neighbor 
spacing. In the next section, we consider the twist-dependent 
changes in cohesive filament contact at the surface of bundles, 
which for sufficiently long bundles evidently counteract the in- 
creases in Eijuik with twist. 



4 Surface energy of twisted bundles 

The ground- state results of our discrete model simulations 
show that while the geometric frustration leads to a cohesive 
energy costs in the bulk of twisted bundles, there is a competing 
tendency for twist to lower this energy, eventually even to val- 
ues lower than the energy of the untwisted state (seen in Fig. [7]). 
In this section, we show that this driving force for twist results 
from an overall decrease in the surface energy of the bundle. 
The surface energy accounts for the deficit of cohesive contact 
at the periphery of bundle, which we analyze in the continuum 
limit of vanishing filament diameter to show that twisting de- 
creases the total number and length of exposed filaments at the 
sides and ends of the bundle. We then compare this to the re- 
sults of our discrete model simulations. Finally, we pair the 
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surface driving force for twist in our continuum model to the 
mechanical cost for bending filaments, to determine the ther- 
modynamic preferences for twist in cohesive bundles in terms 
of the bundle aspect ratio and filament stiffness. 

4.1 Surface geometry of twisted bundles 

In this section, we analyze the surface geometry of twisted fil- 
ament bundles in the continuum limit where filament diameter 
is small compared to both filament length and the lateral size 
of the bundle, with the goal of developing analytical formula 
for the dependence of surface energy on twist. In this spirit, we 
make the additional approximations that: 1) the filament pack- 
ing is locally hexagonally-close packed, with a density that is 
independent of twist; and 2) the shape of the bundle is axisym- 
metric, with an outer cylindrical radius, R. The first approxi- 
mation is clearly violated in the neighborhood of disclinations 
that enter the packing at finite twist. Though, for large bundles 
N I, the local packing is non-hexagonal for only a relative 
minority of filaments. Hence, we assume that the local filament 
spacing and occupied volume fraction change only modestly in 
the bulk of twisted bundles, which is consistent with the density 
of maximally-compact bundles studied in ref. For straight 
filament bundles, the second approximation (cylindrical bun- 
dle symmetry) clearly fails to capture the hexagonal faceting of 
the bundle sides. However, our simulations show that bundles 
become more axisymmetric at high twist as the packing trades 
high-energy corners at the bundle surface for excess disclina- 
tions in the bundle interior 

Based on these assumptions, we now consider the change in 
the surface shape with twist (shown in Fig. |9]). As twist in- 
creases, the helical tilt of filaments away from the center of 
rotation increases as 6{p) = arctan(^lp). Since the contour 
length of filaments are fixed to L, twist requires a change of the 
height, H, the extent of a filament along the pitch axis, accord- 
ing to 

//(p) =Lcose(p) = ^ (18) 

Hence, for non-zero £1, the vertical profile of the ends varies 
with radius. Assuming filaments distribute the taper equally 
over both of the free ends of the bundle, the shape of this ta- 
pered profile is described by H{p)/2, as seen in the example 
bundles in Fig. |9] It is important to point out that the curved 
shape of this profile is unrelated to the geometry of the bundle- 
equivalent surface which encodes the metric properties of in- 



terfilament spacing (see section [33] ). 

Local decreases in the height of the bundle with twist im- 
ply that the lateral radius must necessarily increase in order to 
preserve a constant volume and density. Given and an outer 




Figure 9 Filament bundle radius, R/Rovs twist angle, 6 from 
numerically simulated ground states. The black line is the continuum 
model prediction of eqn \19\ , while the discrete model is represented 
with three bundle sizes: small N = 34 (red), medium N = S2 (blue), 
and large N =1^4 (green). Both the continuum approximation 
(opaque pink surface, radius is + J to account for the filament 
diameter), and the discrete model representations are shown together 
for three example twist angles for a bundle of N = 46. The height of 
the outermost filaments is shown for example C. 



bundle radius, R, the volume within a bundle is easily com- 
puted as V{Q.,R) = 2nQ.-^L{^/l ^ {Q.R)^ - l). Assuming an 
untwisted bundle radius of Rq, the ^1-dependence of R is deter- 
mined from the solution to V{C1,R) 



TtRlL, 



R = Rox l^{^Ro/2)^. 



(19) 



^This follows from the fact that the sum of topological charge interior to 
the bundle and the net deficit of nearest neighbor contacts at the boundary are 



constrained to obey, 2 + L„ (4 — n)Vh (n) = 6. See, e.g 



This formula, while derived from global considerations of vol- 
ume conservation, implicitly encodes the same constraints of 
lateral filament-packing in twisted bundles described by the 
mapping to the dome-like surface: twisting a bundle reduces 
the number of filaments that can be packed at a given radius 
p by cosO(p), implying that filaments must be redistributed 
to larger radii, and notably compares quantitatively with simu- 
lated bundle radii, as shown in Fig.|9] 

4.2 Continuum-limit surface energy 

We now proceed to analyze the contribution to the cohesive en- 
ergy deriving from the twist-induced changes of filament con- 
tact at the surface of the bundle. The surface energy per unit 
area, L, accounts for the loss of favorable cohesive interactions 
due to exposure of non-contacting filament length at the sur- 
face. Simply put, the surface energy attributes a loss of cohe- 
sive energy, 7o/2, per unit length of lost neighbor contact, rel- 
ative to the hexagonal packing of the bulk. Following standard 
arguments for surface energy, the factor of 1 /2 follows from 
the fact that the separation of one contacting pair creates two 
non-contacting filaments 

The surface of a twisted bundle is composed of exposed fila- 
ments that lack the full complement of neighbors to achieve the 
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maximal cohesive energy density. In a twisted filament bundle, 
non-contacting filament lengths arise in two ways. First, fil- 
aments at the radial sides expose lengths of non-contact along 
the entire outer contour. Second, due to the finite contour length 
of filaments, twist leads to "slip" of filament pairs at the ends of 
bundles. The cohesive energy cost in both cases may be derived 
by considering the creation of non-contacting filament length 
by introducing a planar "cut" through a bulk hexagonal array of 



filaments, as shown in Fig. 10 i. The per area exposure of non- 
contacting filament length is determined by the unit normal n to 
the planar cut (the normal of the free surface element) and the 
local orientation of filament tangents, T. Consider, for exam- 
ple, the loss of contact, i.e. the slip length, 4, for the neighbor 
pair shown in Fig. [TO]:, where the tilt of the cut direction is 
along the neighbor separation. In this case, it is straightforward 
to relate the length of surface separating filament ends, ds, to 
slip length, 4 = <i^|sin0|, where | sin0| = |T x n|. In gen- 
eral, summing over the slip of nearest-neighbor contacts yields 
a surface energy per unit area of the form. 



Z = a^|Txn| 



(20) 



where Zq = 7o/^/, and a is a numerical coefficient deriving 
from the orientation of neighbor directions with respect to the 
surface element. For example, for the low-energy sides of the 
bundle where filaments are perpendicular to the exposed sur- 
face (|T X n| = 1), it is straightforward to show that a = 2, due 
to the two fewer neighbors for filaments at the surface relative 
to the bulk. A more detailed calculation (given in Appendix |B]) 
that averages the slip-cost of a cut hexagonal array with respect 
to all cutting directions yields a = Ay/?>/7t c^2.2. For the re- 
mainder of the article, we neglect the the variation in relative 
geometry of neighbor directions and cut directions and simply 
take a = 2 for all bundle surfaces. 

For long filaments, the radial sides of the bundle carry most 
of the surface energy as filaments are normal to the free surface 
along their lengths and |T x n| = 1 is maximal. We define Eside 
as the change in surface energy relative to the untwisted bundle. 



which we calculate using eqns ( 18 ) and ( 19 ) 



^side " 



R/Ro 



(21) 



where Aq = ItiRqL is the side area of the untwisted bun- 
dle. Since the height of bundle side, H{R), decreases more 
rapidly with twist than the lateral growth in radius, Eside is 
a decreasing function of twist. For small twist, Esidei^R ^ 
1) — (3/8)ZoAo(^17?)^, indicating that any finite measure of 
twist reduces surface exposure filaments at the sides. Since the 
length of non-contacting filaments at the boundary is fixed to L, 
this change must derive from a change in the number of surface 
filaments. This demonstrates that a twisted bundle contains a 




Figure 10 (a) A bulk collection of filaments with a cutting surface, 
(b) Zoomed in view. Oblique (c), and side (d) views of the surface 
cut. The orange segments represent the lengths of filaments that are 
now no longer interacting with the neighbors in front of them. This 
length, is, is dependent on the angle, 0, between the surface normal, 
n, and the filament tangent vector, T. This cut corresponds to an end 
surface area section, such as the example shown in (e). 



larger proportion of its filaments in the interior than does an 
untwisted bundle. 

While twist reduces Eside, this reduction comes at the ex- 
pense of increasing surface exposure of non-contacting length 
at the ends of the bundle. We evaluate the surface energy con- 
tribution from one of the ends of the bundle, E^nd, beginning 



with eqn (20). The normal and tangent vectors are defined as 



-H'/2p 



■- cos0(p)z + sin0(p)0. 



(22) 
(23) 



Integration over ends of the surface area of ends, for which 
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dA = 2ndp pi/l + (/fV2)2 (where H' = dpH), yields 
Ee„d = Lo f dA\Txn\ 

J end 



1/2 



1 



, , . , T • (24) 



Analysis of the integrand of eqn ( [241 ) reveals that E^nd has two 
analytically tractable limits whose form depend on the relative 
magnitude of Q,L, which is proportional to the number of he- 
lical turns of a bundle, and sec^ = 1 + {Q,R)^. In the limit 
of infinite length (and finite twist) the surface energy per end 
becomes, 



lim Eend 



7rm"^Zo [arcsinh(tan e) - sin 6] , (25) 



which has a small-^1/? behavior E^ndi^L sec^O) 
^^{QR)^^ In this limit, the surface energy of ends derives 
predominantly from radial slip of neighbor filament pairs ex- 
tending to different heights, for small twist, 4(rad) ^ d\H^\ ~ 
dQ?RL, per pair (see Fig. Hi). In the opposite, limit of vanish- 




Figure 11 An interacting pair of filaments with two possible 
orientation, (a) Radial slip - both filaments share the same azimuthal 
position, (j). Here, the inner filament loses contact at both ends, (b) 
Azimuthal slip - both filaments share the same radius, p. Here, each 
filament experiences reduced contact length at one end of the bundle. 



ing length the surface energy takes the form. 



lim Eend 



7iRQ~^I.o 



cote 



arcsinh(tan0) 



tanO 



(26) 



which has a small-^17? limit, E^ndi^L <^ I) ^\aR\. The 
end surface energy cost in this limit (few helical turns per bun- 
dle) is dominated by interfilament slip between azimuthally- 
separated neighbors, for which 4(azi) ^ \T • ^\d d\Q.R\ at 
small twist (see Fig.pTf)). 

"^Strictly speaking, this is the limit L > Qr^ > Rq, since ^2/? > 1. 



Notably, the ratio of the surface energy contributions cap- 



tured in eqns ( [25) ) and ([26)), which derive from the two distinct 
modes of slip occurring at the ends of twisted bundles, are con- 
sistent with the relative magnitudes of radial vs. azimuthal slip 
in weakly twisted bundles, 4(rad)/4(azi) ~ |0|(L/7?o). This 
implies that the aspect ratio of the bundle, L/Rq, is a key pa- 
rameter governing the twist-dependence of surface energy. In 
Fig. [12] we plot the total surface energy, Esurf = Eside + '^E^nd^ 
as a function of twist angle 6 for aspect ratios ranging from 
L/Ro ^ to the infinite length limit, L/Rq ^. 
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Figure 12 Filament bundle surface energy vs twist angle for various 
aspect ratios of L/Rq. 

First, we consider the behavior of the L/Rq oo limit, 
where the decrease with twist in surface contact at the sides 
of the bundle overwhelms the additional cost of radial slip de- 



scribed by eqn (25 ), such that the surface energy is unstable to 
twist and obtains its minimum for the maximally twisted state 
6 90°. In the opposite limit of L/7?o ^ 1 (vanishing aspect 
ratio), the side and radial-slip contributions at the ends con- 
tribute negligibly to Esurf since both scale with contour length 
of filaments. Thus the surface energy is dominated by pos- 
itive cost of exposing azimuthal slip at the ends of the bun- 
dle, eqn ( [26) ). For small twist, this cost increases linearly as 
Eend {L/Rq <^ l)Mo ~ \^\{^o/L) and diverges for large twist, 
Eend{L/Ro <^l)/Ao ^ coi^/^ e{Ro/L) such that the minimal 
surface energy corresponds to the untwisted state. For interme- 
diate aspect ratios we see a transition between these extremes, 
with the minimal surface energy occurring at = 0° for short 
bundles, and jumping to finite and large values of twist above a 
critical aspect ratio L/7?o — 5.07. 

We assess the quality of the continuum-limit surface energy 
analysis by direct comparison to the numerical simulations of 
the discrete filament model. Simulations of bundle cross sec- 
tions in section 3.1 are carried out in the L/7?o ^ °° limit so that 
surface energy changes with twist are derive only from E^ide 
and the radial-slip contributions to E^nds^ which are both pro- 
portional to L. To extract the surface energy of bundle sides, we 
calculate the excess energy of surface filaments due to fewer fil- 
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ament neighbors than the predominantly six-fold packing in the 
core, 



-^side 



i:r(A.7)- 



-6y{d) 



(27) 



where / G b refers to filaments at the surface of the bundle. Fig. 
13 compares the relative change in surface energy at the sides 
of small, intermediate and large bundles in our discrete model, 
to the continuum expression for Eside is eqn (27) ^ While 




Figure 13 External filament cohesive energy vs twist angle for three 
selected bundle sizes: small N = 3A (red), medium N = S2 (blue), 
and large N = IS4 (green). Black line is the continuum model 
expression eqn pT) . Inset shows number of external filaments. Next 
vs twist angle. 

the circular approximation of the faceted boundary shape for 
straight bundles leads to an underestimation of the surface en- 
ergy change, we find that the continuum expression for Eside 
effectively captures the shape and magnitude of surface energy 
decrease as the bundle is twisted. 

Though not considered explicitly, the discrete model simu- 
lations of the previous section do implicitly count the cohesive 
energy loss due to radial slip at the filament boundaries, deriv- 
ing from the curvature dependence of the cohesive energy. In 
the continuum limit where /Cy J <C 1 we may approximate the 
curvature-dependent prefactor in eqn (p3]) as 



L(l- 



-1/2, 



iL-£s{ij)/2, 



(28) 



where 4(u) — 

DeltUij • Nj = Lij — Lji is the difference in contacting length of / 
with j and contacting length of j with /[^ Note that the implicit 
loss of contact from radial slip deriving from the curvature de- 
pendence does not account for the additional twist-dependent 



^For the discrete model, the value of 6.7O70 was used in place of the 6Y{d) 
in eqn {27) , corresponding to the cohesive energy for nearest and next nearest 
neighbors of a bulk filament in a hexagonal packing of spacing, d. 

^This identity derives from the mapping of curve i to the point of contact 
on j, ^j{si) = ^i{si)+ /^ij{si). Since d/^ij/dsj • = • A/ 

-1 



slip between azimuthally separated pairs (e.g. Fig.[TT])), which 
enters explicitly into Ltj. In the following section we gener- 
alize our discrete model to include these additional costs. We 
derive the surface contribution from loss of filament contact at 
the ends of twisted bundles of infinite length as 



(29) 



where again we take / to be the outer filament of the pair so 



that Lij = L. Fig. 14 compares this surface energy formula- 



tion applied to our small, intermediate, and large discrete model 



bundles, to the continuum model prediction of eqn (25). This 
shows strong agreement over a large range of twist angles as 
the number of filaments per bundles grows sufficiently large. 



As described in Sec. 3.4 we calculate the bulk cohesive energy 




\dRj/dsi\ = (1 + KjAij-N^)) mdLjiC^Lij{l - 



•Nj)). 



Figure 14 End surface energy for three selected bundle sizes: small 
N = 34 (red), medium N = S2 (blue), and large N =IS4 (green), in 
the infinite length limit. Black line is the continuum model 
expression eqn |29j). 

shown in Fig. |8lby subtracting the the surface contributions 
given in eqns p7] ) and ( [29| ) the total, discrete model energy, 

that is, Ehulk = Etot — E sides — '^Egnds- 

4.3 Thermodynamics of surface energy 

The continuum analysis in the previous section demonstrates 
that for sufficiently long bundles (L > S.OIRq) the cohesive en- 
ergy at the surface prefers finite and generically large values of 
bundle twist. Along with the constraints and costs of packing 
frustration at the bundle core, the additional mechanical cost of 
filament bending competes with the surface energy preference 
for twist. In this section, we analyze an albeit simplifed model 
of the twist dependence of cohesive bundles that includes only 
the bending energy and surface energy. Our goal is to deter- 
mine the optimal twist geometry of filament bundles in the ab- 
sence of the additional costs of filament packing at the core that 
were considered in section|3] and which are considered again in 
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the next section. The continuum analysis of these two compet- 
ing energies suggests that minimal-energy bundles generically 
exhibit a degree of spontaneous twist that is highly sensitive to 
both bundle aspect ratio and filament stiffness. 

The mechanical cost to bend a straight filament into a helical 
shape is simply where B is the bending modulus, or 

stiffness, of the filament. In the continuum limit, we compute 
the total bending energy of the filaments in a twisted bundle by 
integrating over the cross- sectional area of the bundle 



^bend 



dN 
dA 



(30) 



where dN /dA is the areal density of filaments in the horizon- 
tal cross section of the bundle. Following volume-conservation 



considerations similar to section |4.1[ assuming a volume frac- 
tion of filaments in the bulk of the bundle, the density at an area 
element located at p is reduced by twist, according to. 



dN 



no 



dA ^l + (^lp)2' 



(31) 



where n^ ^ = {^/3/2)d^ is the cross- sectional area per filament 



in a dense hexagonal packing. Combining both eqns pO\ 
([3T]) to perform our integral in polar coordinates yields. 



^bend ' 



TlBnoL 



(TT(^W 



and 



(32) 



Note that the form of E^end is not an explicit function of fil- 
ament number or Rq, but instead depends only on twist an- 
gle 6 = arctan(^17?). In the limit of weakly twisted bun- 
dles, bending energy exhibits a soft dependence on twist 
limQR^i Eijend = {7iBnoL/4){Q.R)^, while in the limit of large 
twist, the bending cost asymptotically approaches a constant 
value lim^^/^»i Etend = (TiBnoL/S) [2 - 3/{Q,R)] 

We define the total continuum energy as the sum of surface 
and bending energy contributions 



" Esurf ^bend ? 



(33) 



where again Esurf = Eside + '^E^nds^ as defined by the continuum 



expressions eqs. ( [21] ) and ( [24| ). To compare the two types of 
energy, we define a length scale. 



A =^o^/Zo, 



(34) 



which parametrizes the relative costs of bending to cohesive en- 
ergies in filament assemblies. Optimizing Econt with respect to 
the twist angle for fixed Rq, L, and A, we compute the diagram 



of state, shown in Fig. 15 



The phase diagram divides into three principle regions: 
untwisted (L/Rq <C 1), bending energy dominated (^/Rq ^ 
\\L/Rq ^ 1), and cohesive energy dominated {X/Rq ^ 
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Figure 15 Twist diagram of state of for the continuum energy, Econt^ 
showing three regions of optimal twist behavior: Untwisted (left of 
white line), bending energy dominated (upper right "Low Twist" blue 
region), and cohesive energy dominated (lower right "High Twisf'red 
region). The color represents the preferred value of 0, ranging from 
0° (blue) to 90° (red). 



l;L//?o <C 1). I n the untwisted region, twisting a bundle is al- 
ways unfavorable, due to the combined cost of slip at the ends 
of short bundles and bending. In the limit of highly flexible fil- 
aments, where A/7?o <C 1, the transition between untwisted and 
twisted bundles occurs at L/7?o — 5.07; and as filament length 
grows, the balance of side and end surface energy yields an op- 
timal twist that diverges with aspect ratio as Q.R ~ {L/R{)Y^^ . 
For larger stiffness, X/Rq ^ 1, the more significant cost of 
filament bending shifts the boundary between untwisted and 
twisted states to larger aspect ratio sls L/Rq ^ {X/RoY^^. In 
this bending energy dominated region, there is a significant me- 
chanical cost for bending filaments, however, for sufficiently 
large L/Rq, Econt is minimized by a modest twist of a few 
degrees. This arises from the fact that for small twist the 
cost of bending grows as ~ B{Q,R)^, while the surface en- 
ergy decreases as ~ Zo(^l/?)^, leading to an optimal twist of 
{Q.R) ~ {Ro/Xy^^. In this regime, the lowest energy state is 
nearly independent of the bundle length because the bending 
cost dominates the end effects of the surface energy term. Con- 
versely, the amount of twist in the cohesively-dominated region 
(L/RO > l;A/7?o) <C 1) is largely only dependent on L/Rq. 
These two regions are separated by a first order transition for 
bundle lengths beyond a critical size L/Rq > 9.9. In the infi- 
nite length limit, the optimal twist angle jumps from 90° 
to 14.0°, at X/Rq = 2.996. For lower aspect ratios, this first 
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order transition disappears, and the high and low twist energy 
minima merge into one. This critical point is shown as the blue 
dot in Fig.[l5]at L/Rq = 9.9, and A /Rq = 0.425. 

To summarize, we find that the balance of cohesive energy at 
the surface of sufficiently long bundles (L > S.OIRq) and flexi- 
bility favors large bundle twist. For bundles of rigid filaments, 
surface energy drives a more modest degree of spontaneous 
twists in minimal energy bundles above a critical aspect ratio 
that grows with filament stiffness. 



5 Optimal twist of ground-state bun- 
dles: Finite stiffness and length 

The previous sections have identified two competing effects of 
twisted filament bundle geometry: packing frustration of fila- 
ments in the core, and the surface energy cost of non-contacting 
filaments. In this section, we reexamine the energetics of our 
discrete filament bundle simulations, including the full costs 
associated with filament bending and loss of contact in bundles 
of finite length filaments. As described in section [2j cohesive 
energy between filament pairs derives from the local contact 
geometry in twisted bundles. However, for the case of finite 
length filament bundles, the contact length of filament / to 7, 
Lij used in eqn ( [T5] ), must be calculated explicitly to account for 
azimuthal slip at the ends of the bundle (see Fig. [TT])). Defin- 
ing the ends of filaments to be at positions sj ±L/2, Lij is cal- 
culated for any given length in terms of the contact function 
defined in section [2] as 



Lij = Siis'^j = +L/2) - Si{s) = -L/2), 



(35) 



where we follow our original notation that s'j is the arc length 
coordinate of filament j that is the point of contact with Si. It 
can easily be shown that explicitly determining contact length 
for a filament pair along with the curvature-dependent correc- 
tion to cohesive energy, properly accounts for the cost of both 
azimuthal and radial slip of filaments at the bundle ends. 

Evaluating the total energy characterized by an aspect ra- 
tio L/Rq and a finite stiffness corresponding to A/7?o, we de- 
termine the optimal (energy-minimizing) value of twist. The 
phase boundaries separating untwisted and twisted ground 
states for small, intermediate, and large bundles are shown in 
Fig. [16] Importantly, we find that the cost of core packing. 



excluded from the continuum analysis of the previous section 



(see Fig. 15), in combination with the bending cost, elimi- 
nates the regime of stable, weakly-twisted bundles in the bend- 
dominated region, A/7?o ^ 1, that is predicted by in the contin- 
uum model. The discrete model of cohesive filament bundles 
exhibits only two well-defined phases: untwisted bundles for 
sufficiently short or rigid filaments, or highly twisted with a 
twist angle of > 75° for long and flexible filaments. The bun- 



10 



^ 1 



0.1 




Untwisted 




5 10 50 100 500 1000 

L/Ro 

Figure 16 Boundaries of the preferred state of filament bundles for 
three select sizes: small N = 34 (red), medium N = S2 (blue), and 
large N = IS4 (green); overlaid on the continuum model results from 
Fig. [15] Above these boundaries bundles prefer an untwisted state, 
while below, the bundles can lower their energy by adopting a highly 
twisted state (6 > 75°). 



dies in this highly twisted state have the maximum topological 
charge of 2 = 6. 

In the limit of infinitely flexible filaments {X/Rq <C 1), we 
find that the inclusion of core energy shifts the predicted critical 
aspect ratio required for twist from L/Rq > 5.07, to L/Rq > 
8.9. In the other limit of infinitely long filaments (L/7?o ^ 1), 
the critical value of A/7?o, above which the filaments are stiff 
enough to resist the surface energy drive to twist, is X/Rq > 
0.63, A/7?o > 1.41, and A/7?o > 1.78, for = 34, = 82, and 
N = 184, respectively. 

These results show that the additional costs of packing frus- 
tration in the core significantly offsets the gain in cohesive en- 
ergy derived from the surface of long bundles. We find that 
increasing the number of filaments generically increases the 
range over which minimal energy bundles are twisted, substan- 
tially raising the threshold stiffness for the boundary between 
twisted and untwisted states. As a final comment, we note 
the appearance of highly-twisted ground states of our discrete 
simulation model, above the first order line separating highly- 
twisted from weakly-twisted bundles in the oversimplified con- 
tinuum model predicted by optimizing Econt alone, which de- 
rives specifically from the underestimation of Eside (and the 
driving force for twist) in the continuum approximation (see 



Fig. 15). 
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6 Conclusions 

Based on a detailed analysis of discrete numerical model and 
continuum surface model for the cohesive energy in twisted 
bundles, we find that lowest energy state for a bundle of suf- 
ficiently flexible and long filaments is generically twisted.The 
decrease in energy with twist derives primarily from decrease 
surface energy of bundles, which ultimately accounts for the 
fact that twisted bundles expose fewer filaments at the high- 
energy sides than straight bundles. 

While the preference for twist is driven by effects at the 
boundary of the bundle, the complex evolution of cross- 
sectional packing — as evidenced by the universal increase in 
topological charge of the packing — plays a critical thermo- 
dynamic role in stabilizing twisted bundles. Above a criti- 
cal threshold of twist Q 22°, excess 5-fold disclinations are 
needed in the ground-state packing to screen the elastic ef- 
fects of the packing frustration generated by twist. Contin- 
uum elasticity arguments have show that twist decreases inter- 
filament spacing between azimuthally-separated neighbors by 
an amount proportional to {^p)^, ultimately leading to an in- 
crease in energy (per unit volume) that grows as (QR)^ for 
defect- free bundles Hence, in the absence of topological de- 
fects which act to "neutralize" the stresses generated by twist, 
the elastic cost of twisting defect-free bundles grows would 
continue to grow unmitigated at large twist angle, likely over- 
whelming the gains in cohesive energy at the boundary. 

We demonstrate the importance of achieving the appropri- 
ate defect configuration for stabilizing twist by considering a 
class of kinetically-constrained bundle packings in our numer- 
ical simulation model. Unlike our numerical search for ground 
states described in section [3] which explored an ensemble of 
in-plane packings at each value of twist, in Fig. 
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we ana- 
lyze the energetics of simulated packings achieved in the fol- 
lowing kinetically-constrained algorithm. Beginning from the 
energy-minimized packing of an untwisted bundle, we increase 
in small increments of 0,001J. For each Q., we perform a 
steepest-descent minimization of bundle energy (in the L ^ oo 
limit) based on the starting positions of the previous, smaller 
value of twist. Fig. 17 compares the energy of kinetically- 
constrained bundles to the simulated ground states presented 
For small Q.d, the cohesive energies of both 



3.4 



m section 

states are identical as expected since no large-scale filament re- 
organizations are required. This persists until a high enough 
twist forces the ground state bundle into a new configuration at 
Q.d ^0.1. At this point they cannot reach this new ground state 
packing as it requires a global rearrangement of the filaments. 
For example, as shown Figs. |3] and |4j at the transition from 
2 = to 2 = 1 packings the minimal-energy packing changes 
from a defect-free hexagonal packing to packing possessing an 
excess 5 -fold defect near the bundle center. In the kinetically- 
constrained packings, such defects only enter at the boundaries 
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Figure 17 Lower thin line shows the thermodynamically lowest 
energy ground states vs twist rate Q.d (same as shown in Fig.[7]). 
Upper thick lines represent the constrained ground state energies. 
TV = 82 for both. 



of the bundle, migrating only slowly towards the center upon 
further twisting. Fig. |7] shows that further increase in twist 
eventually does allow the constrained bundle to overcome the 
local energy barriers of filament arrangement; however, the to- 
tal energy of states continues to exceed the ground state pack- 
ings inhibiting the stability of the twisted bundle relative to the 
^1 = state. 

The discrepancy between ground state packings and this sim- 
ple model kinetically-constrained bundles demonstrates two 
key points about disclinations in twisted filament bundles. 
First, changes in surface energy will only stabilize twisted bun- 
dles over straight bundles provided that the appropriate, energy- 
minimizing configuration of disclinations punctuates the cross- 
sectional packing. And second, the kinetic pathways by which 
topological defects enter into and migrate throughout the cross- 
section may place strong constraints on whether an externally 
or intrinsically twisting system of filaments is able to achieve 
an optimally-twisted state. This points to the need for further 
work on the mobility of disclinations and plasticity in twisted 
filament arrays. 

Assuming that the assembly kinetics of bundles accommo- 
date the appropriate number and distribution of disclinations, 
we find that regimes of thermodynamically preferred bundle 
twist are separated from untwisted bundles by characteristic 
measures of bundle size. As shown in section |5j our discrete 
model calculations suggest that equilibrium bundles are sponta- 
neously twisted when L > IORq and Ro> ^. Unlike the aspect 
ratio, which is a purely geometric parameter, A varies with in- 
trinsic properties of filament stiffness and cohesive forces. We 
conclude with a brief consideration of the value of A for three 
distinct filamentous systems, with the goal of assessing the 
thermodynamic stability of each to bundle twist: (i) capillary- 
condensed arrays of micro- and nano-fabricated pillars; (ii) car- 
bon nanotube ropes; and (iii) DNA bundles condensed in the 
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presence of polyvalent counterions. We quantify cohesive ten- 
dency for twist in terms of Nc = (A/J)^, roughly the number 
of filaments needed to stabilize surface-driven twist. 

( i) Capillary condensed filaments - On the upper end of fil- 
ament diameters, we find filament arrays held together by cap- 
illary forces, as occurs when filament arrays are drawn from 
a wetting into to non- wetting fluid medium^ ^. In such cases, 
inter-filament cohesion is mediated by liquid bridges spanning 
neighboring filaments in the array. As such, we expect the sur- 
face energy of the bundle, Lq, to be proportional to the sur- 
face tension between wetting and non- wetting fluids, of order 
~ 10 mN/m^ Since stiffness is a strong function of fil- 
ament diameter B ^ d^, bundles of large diameter filaments, 
such as hair^^ and similarly sized glass or polystyrene fibers 
with diameters d > 10 {im and bending stiffness in the range 
B ^3 mN mm^ to 3000 mN mm^, are relatively stable to twist, 
only twisting for filament numbers larger than Nc ~ 3000. On 
the other hand, arrays of more slender nano-fabricated pillars 
d ^ 300 nm are relatively easy to twist by capillary forces even 
for bundles of just a few filaments, as A/^ ~ 1, consistent with 
observations of ref . . 

( ii) Nanotube ropes - On the opposite end of the size spec- 
trum are ropes of single- walled carbon nanotubes, with diam- 
eters typically in the range of J ^ 8 A (for (6,6) SWNTs) and 
d^21 k (for (20,20) SWNTs). Ropes of carbon nanotubes 
are prepared by a variety of methods, from the electric-arc dis- 
charge of graphite to nano-textiles spun from grown nan- 
otube mats^^. Nanotube ropes are typically formed in the limit 
of extreme aspect ratio, L/Rq ^ 100. While the influences 
of tube chirality, metallic/semi-conducting and polydispersity 
properties of nanotubes complicates the simplistic treatment 
of inter-filament cohesion considered here, twisted structures 
have been reported in bundles containing at least tens of the 
SWNT^^^^. Depending on the nanotube diameter estimates 
for stiffness vary considerably, ranging from 5 ^ 30 nN nm^ 
for (6,6) tubes to 4000 nN nm^ for (20,20) tubes, while van 
der Walls attraction between nanotubes in vacuum suggestion a 
cohesive energy per length of 7o ~ 800 pN relatively indepen- 
dent of diameter , from which we estimate A ^ 80 nm and 
2 /im for small and large diameter tubes, respectively. The large 
value oiX/d implies that nanotubes are fairly rigid despite their 
small diameter, presumably due to intrinsic stiffness afforded 
by covalent bonding within tubes. These estimates suggest a 
very modest tendency for nanotube ropes to twist, which varies 
considerably with tube diameter: Nc of 5,000 and 500,000 nan- 
otubes for single tube diameters of 8 A and 27 A respectively. 

( in) Condensed DNA bundles - dsDNA condenses in solu- 
tions of multivalent counterions into tightly packed toroidsSSEll 
and bundles (sometimes referred to in the literature as 
"rods")'SHill. Given a bending rigidity of B ^ 0.24 nN nm^ 
an interaction energy per unit length of 70 6 pN, and 
a center-to-center spacing of condensed DNA, J ^ 3 nm^^ in 



the presence of trivalent cations, we can estimate A ^ 13 nm. 
This sets a critical number of cross-sectional DNA strands to 
stabilize twist as A/^ ~ 18. This result implies that the relative 
flexibility of dsDNa (in comparison to, say, carbon nanotubes) 
is overwhelmed by inter- strand cohesion in nominally sized- 
bundles, and cohesive forces alone may be sufficient to stabilize 
twist in toroidal bundles of dimensions typical for encapsulated 
bacteriophage genomes 50 52 

The simple model estimates above neglect many key aspects 
of inter-filament forces that may further stabilize or inhibit twist 
in cohesive bundles. Notably, the present model does not ac- 
count for the twist dependence of interaction between chiral fil- 
aments, a feature well-known to bias the handedness and drive 
twist of interfilament packing in condensed phases of helical 
molecules from DNA to collagenEEl. Surprisingly, the broad 
range of filament sizes and (achiral) cohesive forces considered 
here suggest that even in the absence of intrinsic or external 
driving torques, thermodynamically preferred twist is the rule 
rather than the exception in cohesive bundles of long and flex- 
ible filaments. 
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A Bundle-equivalent surface 

The geometry of the bundle-equivalent surface can be constructed from simple 
considerations of the space available for packing filaments at a radial distance 
p from the center of a twisted bundle, characterized the length ^(p) of a span 
between points of "self-contact" on a helical curve. In terms of the dual surface 
representation, i{p) corresponds to the circumference of the surface an arc- 
distance p from its pole. This length is determined by considering a helical 
curve at p and the length of the shortest, constant-radius span that connects the 
curve to itself at another point along its distance. The finite length of the span 
between "self-contacts" of a helical filament at p implies that the number of 
finite-diameter filaments that may be packed at a given radius is limited, as is 
the space available for packing finite-diameter discs at a given radius from the 
center of an axi-symmetric surface. Fig. [l^ shows the geometry of this span 
whose length derives from the helix geometry, 

£(p) =Psin0(p) = ^^i==. (36) 



19 



The circumference of the surface grows with p more slowly than 27rp, imply a 
spherically-curved geometry, characterized by a positive Gaussian curvature, 



3Q? 



(37) 



showing that curvature (and hence geometric frustration) is concentrated at the 
pole of that surface (i.e. <^ 1), corresponding to a region near the center of 
twist in a bundle. 




Figure 18 (a) A helical filament (orange tube), (b) The 
bundle-equivalent surface. 

We now consider the relationship between the geometry of the dual surface 
and the topology of the triangular network of nearest-neighbor bonds of discs 
or particles packed on this surface (see e.g. Fig.fTsj)), whose constraints carry 
over to the packing geometry of twisted bundles. For a triangular element, 
connecting three vertices of the packing, the Gauss-Bonnet theorem relates the 
integrated Gaussian curvature of the surface patch within the element to the 
internal angles, Oi, at vertices. 



r ^ 

^patch 1 



(38) 



where we have taken the edges of the patch to be geodesies. Summing this over 
the entire mesh of the packing, we have 



Jmesh „ 



b 



(39) 



where F is the total number of faces in the mesh, V{n) is the number of «-fold 
vertices not at the boundary of the triangulation (each of contributing 2k from 
the sum of interior angles), and 6^ are the interior angles of vertices on the 
boundary of the triangulation. Defining (n) as the number of ^-fold vertices 
at the boundary and using 3F = ^„ [^y(^) -\-{n — l)Vb{n)] , we may rewrite 
eqn {39) as 



dAKG = 27lY,U-l)vin)+Y.[^b 



(40) 



where we note that an ^-fold boundary vertex possesses — 1) interior angles. 
Dividing by 2 tt and making use of the definition of topological charge in eqn 
{16), we have 



where we have defined 

[ dAKo, (42) 

^K J mesh 

which plays the role of the Euler characteristic for a boundary-free surface 
domain, notably increasing as the lateral size of the patch grows large in com- 
parison the curvature radius of the surface (proportional to P). The right-hand 
side of eqn pT) represents distortion of the nearest-neighbor packing from an 
equilateral geometry (0^ / tt/ 3) at the free boundary of cluster, such that the 
deficit between the topological charge of the interior packing and 6x must be 
distributed as boundary distortion of the packing. Approximating the bound- 
ary geometry of the packing as circular and integrating the Gaussian curvature 
within a packing of arc -radius R we find. 



1 



[1 + (n/?)2]3/2 ' 



(43) 



which increases from at small twists as 3{Q.R)^ /2, to a maximum of 1 in the 
limit \Q.R\ oo. By defining Qi^ = as the ideal discHnation charge that 

perfectly neutralizes the integrated Gaussian curvature in eqn ( [4T) , we arrive 
at eqn fTT) in the main text. Specifically, when the actual topological charge 
can achieve the ideal value (Q = Qi^), the distortion of the packing at the outer 
boundary from equilateral geometry (6^ = ^/3) can vanish. 



B Lattice orientational dependence of 
surface energy 

As described by eqn {20) , the energy per unit of exposed surface area of a bun- 
dle, £, is dependent on the angle, 0, between the local filament tangent, T, and 
the cutting plane normal vector, n. This definition separates the surface corre- 
sponding only to the loss in contact lengths between filament pairs, as opposed 
to surface area associated with filament ends. However, as noted in section 
|4.2) , there is an additional dependence of L on the orientation of the cutting 
plane with respect to the lattice directions at the surface, which is defined in 
terms of the component of n that lies in the horizontal cross section, (shown in 
Fig.[l9) 

n-T(n T) 



l-(n.T): 



(44) 




10 



cut by a plane with 



Figure 19 Top view of bulk filaments in Fig. 
the normal, n. Filament / and its neighbor j are specified. The 
interfilament spacing is d, and ^ is the angle between the horizontal 
cross section component of the normal, n^, and the vector Xj — Xj. 



We consider the loss of contact for a single filament, i, for each its six nearest 
neighbors, j, 

4(;) = |tan0|J|cosv/,-|, (45) 

where y/j is the angle between Aij and . Defining the angle of this orientation 
as ^, the smallest angle between the bond directions in the hexagonal lattice 
and nx we have Yj = ^ + ^j/3. Using the fact that area per filament at the 
cutting surface is n^^ sec0 and defining the dimensionless parameter a as in 
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eqn {20) , we have 

«W = ^i:|cos(^ + 7rj/3)|, (46) 
^ j=0 

where tiq — 2/y/3d'^. Summing over the filament neighbors we have, 

a(^) = ^ cos^, for -71/6 71/6. (47) 

v3 

Two limiting cases of a: 1) the lattice vector is aligned with (high surface 
energy), yielding a maximum = 0) = 4/^3; and 2) the lattice vector is 
^ ^ ±k/6 maximally offset from (low surface energy), yielding a mini- 
mum = 7i/6) ^ 2. For surface elements at the ends of the bundle, the 
distribution of ^ roughly visits all orientation equally, hence, suggesting the 
appropriate value of a is average with respect to ^: {a) = 4y/3/K ~ 2.2. 
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